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CLASSICAL FOUNDATIONS OF MANY-PARTICLE QUANTUM 

CHAOS 

BORIS GUTKIN, VLADIMIR OSIPOV 


Abstract. In the framework of semiclassical theory the universal properties of 
quantum systems with classically chaotic dynamics can be accounted for through 
correlations between partner periodic orbits with small action differences. So far, 
however, the scope of this approach has been mainly limited to systems of a few 
particles with low-dimensional phase spaces. In the present work we consider N- 
particle chaotic systems with local homogeneous interactions, where N is not nec¬ 
essarily small. Based on a model of coupled cat maps we demonstrate emergence 
of a new mechanism for correlation between periodic orbit actions. In particu¬ 
lar, we show the existence of partner orbits which are specific to many-particle 
systems. For a sufficiently large N these new partners dominate the spectrum of 
correlating periodic orbits and seem to be necessary for construction of a consistent 
many-particle semiclassical theory. 


1. Motivation and Goals 

“Sadly, searching for periodic orbits will never become as popular as a 
week on Cote d’Azur, or publishing yet another log-log plot in Phys. 
Rev. Letters.” 


— P. Cvitanovic, et al, p] 

Already at the dawn of the quantum era it was realized that properties of quan¬ 
tum systems crucially depend on their classical dynamics. While the eigenvalues 
of integrable systems can be explicitly related to the set of integer numbers by the 
Bohr-Sommerfeld quantisation rules no such regular structure exists for systems 
with complex dynamics. At first sight, the energy spectrum of complex systems 
like nuclei resembles a structureless set of random numbers very much dependent 
on particularities of the system. This, however, turns out to be not entirely true. 
The famous Wigner-Dyson-Mehta conjecture asserts that the spectrum of complex 
quantum systems on the scales of the mean level spacing between eigenvalues is 
universal and can be described by Random Matrix Ensembles within the same sym¬ 
metry class. In particular, all n-point correlation functions of eigenvalues can be 
derived analytically by using Random Matrix Theory (RMT). Indeed, such spec¬ 
tral statistics have been observed in many real and numerical experiments with 
various systems ranging from compound nuclei to complex molecules and atoms. 
Furthermore, it was realized in the early 1980s that even single-particle systems 
generically exhibit the same universality if their classical dynamics are chaotic. In 
the last three decades a substantial progress has been achieved in understanding of 
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the origins of this universality through the application of the semiclassical theory 
to quantum chaotic systems. In particular, the Gutzwillcr’s trace formula allows to 
express eigenvalue density function through the sum of unstable classical periodic 
orbits. By using this representation the correlations between system eigenenergies 
can be straightforwardly related to the correlations between actions of periodic or¬ 
bits. This approach was pioneered by M. Berry in the seminal paper [3], where the 
diagonal correlations between periodic orbits were taken into account. It was also 
noted there and in subsequent works that in order to obtain the full RMT 

result one would need to include non-trivial correlations between periodic orbits as 
well. 

On the quantitative level the non-trivial correlations between periodic orbits were 
first taken into account in the groundbreaking work of M. Sieber and K. Richter [7]. 
On the basis of classical chaotic dynamics they demonstrated existence of periodic 
orbit pairs with close actions. Generally speaking, such orbits traverse approxi¬ 
mately the same points of the configuration space but in a different time order. By 
including the Sieber-Richter pairs, and their natural generalizations it turned out to 
be possible to derive the full RMT result for universal spectral correlations [6]]. 

So far, however, this remarkable progress has by and large been restricted to the 
systems composed of just few particles. Although, formally Hamiltonian system of 
N particles in d dimensions can be thought as just one particle in Nd dimension, such 
dynamical interpretation does not allow to use automatically the “single-particle” 
semiclassical theory when N grows simultaneously with hr 1 . This is because the 
standard semiclassical limit assumes a fixed dimensionality of the system while the 
effective Planck’s constant h e ff tends to zero. Another, often considered “semi- 
classical” theory corresponds to the thermodynamic limit, where N —>■ oo but h e ff 
is fixed. As can be clearly seen on the example of bosonic systems, this limit es¬ 
sentially differs from the one where N is fixed and h e ff — > 0 |I0J [IT]. Indeed for 
interacting bosons the two limits correspond to very different classical Hamiltonians 
associated with the first and the second quantisation of the problem, respectively. 
At the present, very little is known about the limit where the number of particles N 
and the effective dimension of the one-particle Hilbert space grow simultane¬ 
ously. This is especially frustrating considering that the original Wigner’s conjecture 
relates to many particle systems. One might wonder whether the “single-particle” 
correlation mechanism between periodic orbits is still relevant for A r -particle chaotic 
systems. Do these correlations lead to a spectral universality, when N is large? At 
which energy scales? Does the answer depend on the details of the Hamiltonian? 

The central goal of the present paper is to show that A^-particle chaotic systems 
with local nearest-neighbor interactions allow a more general correlation mechanism 
between periodic orbit actions then one which is known for one particle systems. By 
using the one-particle interpretation of classical dynamics it is easy to demonstrate 
that standard partner orbits (e.g., Sieber-Richter pairs) exist for an arbitrary N. As 
expected, these pairs of periodic orbits closely approach each other in the large Nd 
dimensional configuration space. It turns out, however, that completely different 
classes of periodic orbit pairs with close actions exist, as well. In particular, these 
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non-standard pairs traverse completely different points of the large configuration 
space. As we show in the body of the paper, they can be interpreted as two- 
dimensional extensions of standard Sieber-Richter pairs, where periodic orbits sweep 
surfaces rather then lines in the small configuration space of the corresponding 
single-particle system. From the semiclassical point of view the contributions from 
new families of correlating orbits should become dominant over the standard one 
whenever the number of particles N exceeds certain threshold value which depends 
logarithmically on ^T//- We believe that in this regime it is essential to take them 
into account in order to construct a proper semiclassical theory. 

The paper is structured in the following order. In the next section we present 
our main ideas on the heuristic level. The emerging mechanism of correlations 
between periodic orbits is then discussed in the framework of the particular model 
of coupled cat maps. This model is introduced in Section 3 where we also describe its 
main dynamical properties. In Section 4 a special duality relation between periodic 
orbits is introduced. Implications for periodic orbit correlations and few remarkable 
mathematical identities stemming from this duality are discussed here. In Section 5 
we describe general mechanism of correlations between periodic orbits and provide 
exhausting numerical evidence of its existence. We explicitly construct a number 
of periodic orbit pairs with small action differences. A general formula for action 
differences within such pairs of orbits is derived here. The implications of symmetries 
on the orbit correlations are studied in Section 6. In Section 7 we extend our 
results to the perturbed maps. Finally, in Section 8 we discuss the relevance of our 
results for generic many-particle systems with local homogeneous interactions and 
its implications for a semiclassical quantum theory. 


2. Main Ideas 


In the present work our attention is focused on chaotic systems of W-particles 
with local homogeneous interactions. To be more specific, let us write a prototype 
Hamiltonian which governs the dynamics of suchlike systems: 

N 2 

(2-1) T~L ^ ' X I" Vi n (lC n mod N j %n +1 mod TV") T V (x n ) . 

2m 

n= 1 


Here x n is the particle coordinate belonging to a target space X e.g., X = R d ,S d 
and p n is the corresponding momentum. Since the energy of the system is preserved, 
the dynamics of < \2.1 ) is restricted to the 2dN — 1 energy shell f 2e, determined by the 
conditions T-L(x i.. .xn,Pi ■ ■ ■ Pn ) = E. Note that R is manifestly invariant under 
the cyclic shift 


( 2 . 2 ) 


C\ : n —> n + 1 mod N 


of the particle numbers. In addition to the shift invariance we will assume that the 
dynamics induced by R on Qe are fully hyperbolic. The last requirement implies 
that the derivative of the Hamiltonian flow has exactly dN — 1 positive Lyapunov 
exponents in almost every point of his. Informally speaking, the above Hamiltonian 
describes a closed string of particles uniformly interacting with their neighbors and 
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FIGURE 1. In the upper figure (a) is shown the model of Boltzmann gas 
with the Hamiltonian of the form l2.l1 Each ball of the radius a is constraint 
to a particular cell of the chain and interacts with the two neighbors through 
the openings of the size L < a. On the bottom figure (b) is shown the 
model of coupled cat map. Each particle lives on a circle of unit length 
and linearly interacts with the two neighbors at the discrete moments of 
time. The corresponding time evolution of its coordinate q, is described by 
eq. 


(3.4). 


an external potential V such that the resulting dynamics are chaotic. One possible 
realization of such Hamiltonian is provided by a Boltzmann gas type model shown 
in fig. [3}i. 


Remark 2.1. The Hamiltonians of the type (2.1) also appear naturally in connection 


with the spin chain models and interacting bosons on a circular lattice. For instance, 
the Bose-Hubbard Hamiltonian written in the second quantisation language: 


N 


n = 


u 


N 


i— 1 2—1 


UATli 


i) 


N 

i= 1 


rii 


hi = b\bi 


can be cast into the form similar to (2.1) after substitution bi = 4^(fj + ipf), where 


Xi and pi should be interpreted as the coordinate and the momentum operators of the 
i-th harmonic oscillator. 

2.1. Time evolution. Given some initial conditions the time evolution of a N- 
particle system can be represented as collection 

r ( r ) = (7i 72 (t), • • •, 7iv(t)}, t G R 

of one particle trajectories 7 n (r) = (x n (r),p n (r)) G V, n — 1,... N in the one- 
particle 2d-dimensional phase space V := T*X . In what follows we will refer to 
T(r) as many-particle orbit and as many-particle periodic orbit (MPO) if T(r) is 
closed i.e., 7 n (r) = 7 n {r + To), n — 1,... N after a period r 0 . Rather then treating 
evolution as a continues flow, it is convenient to discretize dynamics by introducing 
a discrete Poincare map acting on a section of the phase space. To this end we fix a 
Poincare section V which is a 2 dN — 1 dimensional hypersurface invariant under the 


shift (T 1 , see eq. ( 2 . 2 ), and observe the system at the discrete moments of time r t , 
f 6 Z when many particle orbit r(r) pierces V. For the system shown in fig. 0 ' one 
possible choice for the Poincare section would be to set it in accordance with the 
pairwise collisions of the balls. In this case r t would correspond to the moments of 
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time when balls collide with each other. The above procedure defines the Poincare 
map (Ev acting on the 2dN — 2 dimensional phase space V N = V D VLe, obtained by 
the intersection of V with the the energy shell. 

Under the action of <fr N any MPO T of the system (2.1) can be naturally seen 
in two different ways. First, Y can be thought as collection of the points in the 
IV-particle phase space V N registered at the discretized moments of time: 

5 1 (r) = {r(T 1 ),r(r 2 )...r(T T )}. 


In other words the evolution of the whole system can be seen as propagation of one 
particle in the “large” phase space V N . Accordingly 5 1 (r) is naturally interpreted 
as a discretized one dimensional line. Such a reduction to one particle system is a 
widely used tool in dynamical systems. For example it allows to treat the Boltzmann 
gas as a many-dimensional billiard [T 2 ] . 

The second conceptually different way is to look at Y as a set of N ■ T points in 
the “small” single-particle phase space V: 

£ 2 ( r ) = (7n.il n — 1,... ,N, t — 1,... ,T}, 

where 7 ^ := 7 n {jt) marks the position of the n’th particle at discrete moment of 
times t — 1,... ,T. For the sake of convenience of presentation it is instructive to 
introduce the map T r : Z 2 T —> V which attaches to each point (n,t) of the discrete 
“particle-time” space Z 2 T := {1,..., N} x {1,..., T} the corresponding point 7 ^ 
of the phase space V. With such notation we have 

S 2 (r) = *r(. z 2 t ). 

In a (properly scaled) continuous limit N, T —y 00 the time evolution of the above 
many-particle system can be seen on an intuitive level as the propagation of a one 
dimensional closed string y(£, r) = 7 (£ + £q, t), where £ G [ 0 , £q] measures the length 
along it and r is continues time parameter. Accordingly, any periodic orbit of a 
period r 0 can be represented by a two dimensional torus swept by y(£, r) in V 
during the time r € [0, 7 "o]. This contrasts with the single-particle interpretation 
where periodic orbits of the system are naturally represented by one dimensional 
lines. It should be emphasized that we actually neither require the interaction 
between particles be attractive nor we consider the continues limit. So the resulting 
time evolution of the system might be very different from a string-like motion where 
all particles keep a linear order in the target space. Nevertheless, regarding the 
correlation mechanism between periodic orbits it turns out to be quite useful to 
think about a MPO T as of discretized two dimensional surface 5 2 (T) in the one- 
particle phase space V (rather then one dimensional line 5 1 (r) in the many-particle 
space V N ). 


2.2. One particle partner orbits. Before focusing on the action correlations be¬ 
tween periodic orbits in many-particle systems, let us briefly recall the correlation 
mechanism in one particle Hamiltonians 0 Ej. A sufficiently long periodic orbit 
typically has a number of self-encounters. These are the stretches of the trajec¬ 
tory where it closely approaches itself in the configuration space, see fig. [2] For a 
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trajectory possessing sufficiently long encounters the hyperbolic nature of the dy¬ 
namics guarantees existence of several (the exact number depends on the encounter 
structure) of partner periodic orbits which traverse approximately the same points 
of the phase space, but make different switches at the encounters. The action dif¬ 
ferences between the partners are accumulated primarily at the encounter stretches 
and decrease exponentially with their lengths. Accordingly, the periodic orbits of a 
hyperbolic system can be organized into families, where all members have approxi¬ 
mately the same actions. 




FIGURE 2. On the left is shown the caricature of a Sieber-Richter pair 
with one encounter shown as a (red) box. Such partner orbits exist only 
in systems with time reversal symmetry. On the right it is shown the 
caricature of a periodic orbit (solid line) and its partner (dashed line) with 
two encounters e\, e 2 ■ Their symbolic representations are given by (2.3). 



FIGURE 3. a) Schematic depiction of ID symbolic representation of 
the periodic orbit on fig. UK Two encounter regions correspond to the 
identical sequences of symbols labeled by a and b letters respectively b) 
The structural diagram of the same periodic orbit. The encounter regions 
correspond to two pairs of points on the circle T 1 . After “gluing” the circle 
at these points we obtain the graph shown in the second figure. This graph 
can be also interpreted as Feynman diagram for a zero-dimensional sigma 
model which describes universal spectral correlations in quantum systems 
with classically chaotic dynamics, see [6]. 


2.2.1. Symbolic dynamics. The intuitive description of the above correlation mecha¬ 
nism between periodic orbits can be formalized by making use of symbolic dynamics 
mm- Within this approach the first step is to fix a Poincare section and consider 
the discretized time evolution. Assuming that the reduced phase space of the sys¬ 
tem allows a finite Markov partition |J ^ V a , any trajectory 7 can be encoded by 
a doubly infinite sequence of symbols from some alphabet A of a finite size: 


. . . ( 2 _ 2 &— iCLq.CLiCL2 • • • 
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Here each symbol a t registers the position of the particle with respect to the partition 
at the discrete time t, i.e., a t equals to a G A if 7 (t) G V a . Correspondingly, a n- 
periodic orbit 7 can be associated with a cyclic sequence of the length n\ 

a 7 = [cqGQ ... a n ], a,i G A. 

On the level of symbolic dynamics each l-encounter of the length p in 7 is a string 
of p symbols e = eie 2 .. . e p which appears in a 7 a number l > 1 of times, see fig. [3^. 
Given a periodic sequence a 7 its pair sequence a 7 is constructed by rearranging the 
encounters and stretches of symbols connecting them in a way that any subsequence 
of consecutive p symbols existing in a 7 reappears in a 7 exactly the same number of 
times and vice versa. Since the approximate position of a point in the phase space 
is determined by a finite string of symbols, the above property guarantees that the 
corresponding orbits traverse approximately the same points of the phase space. For 
instance, the following two sequences: 

(2.3) a 7 = [aeide2cei6e2], a 7 = [aei&e2ceide2], 

with the encounters e\, €2 of some length p and stretches a, b, c, d of an arbitrary 
length connecting them provide the symbolic representation of partner orbits shown 
on fig. [ 2 ] 3 . 

Structurally different families of partner orbits { 7 } can be distinguished by their 
diagrams Jy accounting for different order of encounters. Each of the diagrams is 
obtained by substituting encounters of 7 with the points and placing them (in the 
same linear order) on a circle T 1 . The points belonging to the same encounters 
are then identified (i.e., T 1 is “glued” to itself at these points), see fig. & The 
resulting graph Jy carries all structural information on the correlations between the 
partner orbits { 7 }. In particular, it defines the number of partner orbits within the 
family. It worth mentioning that in the semiclassical theory based on the correlations 
between periodic orbits such graphs can be associated with Feynman diagrams of 
the corresponding sigma model [ 6 ]. 


2.3. Many particle partner orbits. We turn now our attention to the iV-particle 
systems with the local type of interactions ( 2 . 1 ). At the beginning of this section 
we argued that in the limit of continuous N the periodic orbits of such a system 
can be regarded as two-dimensional surfaces in the phase space. Taking this as a 
guiding principle let us first analyze on an intuitive level its implications on the 
correlation mechanism between periodic orbit actions. By the analogy with the 
single-particle case, we demand that the two surfaces corresponding to some pair 
of partner periodic orbits T, T pass through approximately the same points of the 
phase space i.e., 

S 2 (T) = T r (Z^ T ) » (Zl T ) = 5 2 (f). 


For this to happen each periodic orbit must have at least one /-encounter, which is 
an element of the configuration space where the surface approaches itself sufficiently 
closely for a number / > 2 times, see fig. [IJ The existence of the encounter implies 
that there are l disjoint domains C Z 2 T ,i = 1,...,/ related to each other by 
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the translation shift ^ : (n, i) —y (n + n! mod N,t + t' mod T): 

(2.4) E® = ■ E w , for i = 1,..., l, n E ij) = 0, for i ± j 

such that all domains are mapped into approximately the same points of the target 
space V N : 

T r (F (1) ) « T r (F (2) ) « • • • « T r (F W ) ~ M(A (1) ) « ^r(£ (2) ) ~ ■ ■ ■ « T f (M). 

Conversely, whenever a closed surface possesses encounters it might be possible to 
find a pair surface traversing almost the same points of the phase space. Therefore, 
the problem of classifying structurally different families of correlating MPOs can be 
on a geometrical level understood as one of finding all topologically non-equivalent 
ways to “glue” a two-dimensional surface to itself. This is very much similar to 
the single-particle case where we need to classify the non-trivial ways to “glue” 
one-dimensional lines. 




FIGURE 4. On the left is a geometric representation of MPO with an 
encounter (shown in red). On the right is shown symbolic representation of 
suchlike MPO. Here the alphabet is composed of four symbols {1, 2, 3, 4}. 
The encounter region is encoded by 2D sequence of symbols (shown in 
green) which reappears in two different places of Z^ T . (color on line) 


2.3.1. 2D symbolic representation. The above intuitive geometrical picture of part¬ 
ner MPOs can be made precise through introduction of two dimensional (2D) sym¬ 
bolic dynamics. Within this approach, it is assumed that each many particle (peri¬ 
odic) orbit T can be encoded by a two dimensional (periodic) lattice of symbols a Uit , 
(■ n , t) G Z 2 , where symbols a n ^ t belong to some alphabet A of a small size. Naturally, 
each MPO T is represented then by N x T toroidal array of symbols: 

A r = {«n,t| CM) e Z 2 t }. 

Remark 2.2. One formal way to construct such a symbolic representation, would 
be to use some partition U«eA °f ^ e one-particle phase space V . At a time t 
a symbol a U)t G A is attached to the n 'th particle, if the point belongs to the 
corresponding part of the partition: V (ln t . A symbolic representation of this type has 
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FIGURE 5. a) Schematic depiction of 2D symbolic representation for a 
MPO r with two encounters. The first 3-encounter has winding numbers 
(0, 0). It is represented by three identical sets of symbols at the frame-like 
domains E^\i = 1, 2, 3 shifted with respect to each other on the torus Z^ T . 
The second 2-encounter is represented by two identical sets of symbols at 
the strip-like domains E^\i = 1,2 and has winding numbers (0,1). b) The 
structural diagram Ap of the MPO on the left figure. The two encounter 
regions correspond here to 3 identical (blue) oval lines and 2 (red) lines 
winding around the torus. The two-dimensional surface J~r is obtained by 
“gluing” the torus along the identical lines. Compare with the one particle 
structural diagram shown in fig. [3| 


been constructed for a model of coupled lattice maps in [22J 

To facilitate further discussion we introduce the map ipr ■ ^nt ~^ which assigns 

to each point (n,t) E the corresponding symbolic representation a n j of T i.e., 

ifa(n,t) = a n>t . 

In what follows the array Ap = -0r(Z^ T ) is referred to as 2D symbolic representation 
of T. By the analogy with one-dimensional symbolic dynamics of hyperbolic systems 
we require few basic properties for such 2D encoding. 

• Small alphabet: The number of symbols in A is independent of N. 

• Locality: The symbols in a neighborhood of each point (n,t) define an ap¬ 
proximate position in the phase space of the n'th particle at the time t. More 
precisely, let p be an odd integer and let 

= {(i, j)|i = n — (p— l)/2 mod N,... ,n + (p — l)/2 mod N, 
j — t — (p — l)/2 mod T,... ,t + (p ~ l)/2 mod T} 

be a square-like set of p x p points around (n,t)- We require that for any 
two MPOs T = { 7 n,i|(n, t) E Z^ T }, T 7 = { 7 ^| (n,t) £ ^ T } having the same 

symbolic representation at the set 7 Zp l ' t ' 1 i.e., ipr= ifr the 
distance between the corresponding points in the phase space is bounded in 
an exponential way: 

\~in,t - i„f \ < V, 


(2.5) 


where A is some constant. 







MANY-PARTICLE QUANTUM CHAOS 


10 


• Uniqueness : For each MPO there is a unique 2D symbolic representation. 
In other words if r(Z^ T ) = ^r'(^ T )i then necessarily T = T'. Note that 
we do not require the opposite i.e., there might be 2D sequences of symbols 
which do not correspond to any real MPO. 

Whether or not there exists symbolic representation of periodic orbits with the 
above properties is a highly non-trivial question, which should be, in principle, ex¬ 
plored for each system individually. It is known that such 2D symbolic representation 
can be constructed for some coupled lattice maps [22j. In the body of the paper 
we also demonstrate the existence of suchlike symbolic dynamics for the model of 
coupled cat maps. Leaving for a while the problem of existence aside let us describe 
how the 2D symbolic representation can be utilized in order to find partner MPOs. 

2.3.2. Partner orbits. First recall that by definition partner MPOs T, F pass through 
approximately the same points of the phase space. The locality property of the 
symbolic dynamics implies that the symbolic representation Ar in a neighborhood 
of a point (n, t ) must coincide with one of its partner Af in the neighborhood of 
some another point ( n , t). This motivates the following definition. 

Definition 2.3. Let p be a fixed positive integer. We call two symbolic representa¬ 
tions Ar, Af as p-close iff any p x p square of symbols: 




/«1,1 

OL 1,2 • 


(2.6) 

[ot]p — 

£*2,1 

Ot 2,2 • 

£*2,p 



K a p,i 

a P) 2 ■ 

£*p,p/ 


Oii : j G M 


appears the same number of times (which might be also zero) in both Ap and Af. 


Clearly if two MPOs T,F have p-close 2D symbolic representations Ap,Ap then 
T and F are partner orbits, with the parameter p controlling how close the sets 
S 2 (T), S 2 (T) approach each other in the phase space. 

Given a 2D symbolic representation Ap a p-close 2D sequence might only exist 
if Ap contains at least one /-encounter of a “width” p. In other words, if T h as a 
partner orbit F there exist l disjoint domains Fa*) c Z 2 t , i — 1,... ,1 satisfying (2.4) 
such that 2D subsequences of symbols at all coincide (see fig. [i]): 

(2.7) ME (1) ) = ME {2) ) = ■ = ME (l) ). 

A general encounter region can be defined by fixing closed (translationally related) 
discrete line contours T E a),i = 1 ,...,/ and expanding them into “strips” of finite 


width. The resulting sets A), i = 1 ,...,/ must fulfill the eq. (2.4) and satisfy the 


following two requirements: 1) The interior of each T E « should contain at least one 
point which does not belong to i.e., E® has to be multi-connected. 2) For 
any (n ,t ) G T E (i) one has TZp' J> C F/b) . Structurally different encounters can be 
distinguished by their winding numbers w = ( W\,W 2 ) specifying how many times 
the path T E a) goes around the torus Z| T in the particle and the time directions, 
respectively. 






MANY-PARTICLE QUANTUM CHAOS 


11 


Let r be a MPO and let Ar be its symbolic representation. In general T might 
contain a number A f of encounters with different multiplicities l n and winding num¬ 
bers w n , n — 1,.. .A f. For each encounter there exist l n corresponding sets of points 
En' } c Z^ T ,i = 1 related by shift translation such that symbolic represen¬ 

tation of T on all these sets coincide, see fig. UK To distinguish between different 
classes of MPOs let us introduce for each T the corresponding diagram E p keeping 
the topological information about its encounters structure. To this end substitute 
Z^ T with a 2 -dimensional torus T 2 and each set En 1 with a line Tn' 1 c T 2 (whose 
exact shape of no importance) having the same winding numbers and the same ge¬ 
ometrical order as the corresponding encounters. We then identify the points of T 2 
which belong to the lines Tn ' 1 , i = 1,... l n , see fig. [HJd. (Informally speaking, we glue 
the torus to itself along the lines from the same encounter). The resulting surface 
Ey will be referred to as structural diagram of T. It has the the same role as the 
structural graph for one-particle periodic orbits, see fig. |3j In particular, Ey defines 
the number of T’s partners. 

The symbolic representation Ap of a partner orbit f of T can be recovered from Ap 
by exchanging the patches of symbols outside the encounter regions. The exact rules 
for such exchange and the number of partner orbits obtained in this way depend 
only on the structural diagram Ey of the MPO. Below we illustrate this procedure 
for periodic orbits with the minimal possible number (one and two) of 2 -encounters 
having winding numbers: w — (0, 1 ), w — ( 1 , 0 ) and w = ( 0 , 0 ). As will be argued in 
Section | 8 j MPOs with such encounters are the most relevant from the semiclassical 
point of view. 
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FIGURE 6. On the left figure (a) are shown two partner MPOs with 
two encounters winding around the torus Zj( T in the particle direction. 
The 2D symbolic representation of the second MPO is obtained from the 
first one by exchanging symbols in the domains B and D. In one particle 
interpretation this is precisely the pair of trajectories shown in fig. §>• On 
the right figure (b) are shown two partner orbits with encounters winding 
around the torus in the time direction. Their construction is precisely the 
same as on the left figure with the switched particle and time directions. 
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FIGURE 7. The pair of MPOs with one encounter having zero winding 
numbers. The 2D symbolic representation of the second MPO (right) is 
obtained from the first one (left) by exchanging symbols in the domains A 
and B. 


a) “Horizontal” encounters: w = (1,0). Encounters of this type wind up torus 
in the particle direction. A simple example is provided by the two horizontal strips: 

E = {fi,... fi + p} x {1,..., N}, E' = {t 2 , . ■ ■ t 2 + p} x {1,..., A}, 

\t 2 — ti\ >j)GN, with the symbolic sequences coinciding on these sets: 

e = me) = Me'). 


If T has two such encounters with the symbolic representations Ei and E 2 , respec¬ 
tively the symbolic representation of its partner f can be constructed simply by 
using the same procedure as in the single-particle case, see eq. (2.3) and fig. [ 6 ^ 1 . 
b) “Vertical” encounters: w = (0,1). These encounters wind up around torus 
in the time direction. For instance the following pair of vertical strips 


E = {1,... ,T} x {ni, • • .ni +p}, E 1 — {1,... ,T} x (n 2 , ...n 2 + An}, 

1 77,2 — 711 1 >p6l caring the same symbolic representation defines an encounter with 
the required winding numbers. The construction of partner orbits here is completely 
analogous to the previous case up to the switch between particle and time directions, 
see fig. i>- 

c) “Annular” encounters: w = (0,0). For the construction of partner orbits 
in this case, it is sufficient to have a periodic orbit with just one encounter. Let 
A and B be two sets of the points from corresponding to the interiors of the 
encounter sets E, E' , as shown in fig. [7j If Ar is the symbolic representation of MPO 
r, the symbolic representation Af of its partner orbit T is obtained by exchanging 
the symbols A = ip r(A) with B = M B ), see fig.0 


2.4. Comparison between one-particle and A-particle partner orbits. Re¬ 
call that a MPO might have a partner orbit only if it poses an encounter. What does 
the existence of an encounter in A-particle systems implies on the level of individ¬ 
ual particles? Regarding a MPO as single-particle closed orbit in many-dimensional 
phase space the standard definition of encounter would imply that all particles si¬ 
multaneously reaper at (approximately) the same positions at least twice during the 
period of motion. This is indeed the correct interpretation for horizontal type of 
encounters having winding numbers (1,0). However, encounters with other winding 
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numbers have a different interpretation in the many particle systems (2.1). For in¬ 


stance, if the two encounter domains have the shapes of horizontal strips, see fig. 
only a number of the particles ni, n\ + 1,..., n\ +p perform approximately the same 
motion in V as the particles u- 2 , + 1,..., 712 + p, but they do it for all times be¬ 

tween 0 and T. Furthermore, having an encounter of an annular type shown in fig. [7] 
implies that a group of particles retrace positions of another group of particles in a 
different point of time. We would like to emphasize that only the horizontal strips 
which wrap the torus in the particle direction would qualify as encounters in the one 
particle interpretation. So the 2D representation of MPOs is absolutely necessary 
in order to account for partner orbits with other types of encounters. 


So far the whole discussion of the correlation mechanism between MPOs has been 
restricted to a pure conceptual level. In the next section we introduce a particular 
model of the type (2.1). Within this model we are able to demonstrate existence of 


the partner MPOs and estimate the differences between their actions. 

3. Model of coupled cat maps 

From the technical point of view it is often convenient to consider dynamical 
systems whose time evolution is given from the start by a discrete map rather 
then by a continues flow. One of the best studied and understood examples of 
systems with chaotic dynamics are provided by cat maps which are the hyperbolic 
automorphisms of the unit 2-torus V = T 2 [20, 16] . To introduce a many particle 
set up we couple N cat maps in a linear way, such that the resulting 2 N dimensional 
map: 

T v ; V N y V N , 


V N = T x T ■■■ x T 

"-V-' 

N 


is again a hyperbolic automorphism of the unit 2N-torus V N = T 
Specifically, the generating function of the map <1?^ is defined as 


2 N 


see 


fig- HI 


N 


N 


(3.1) S(q t , qt+l) d ^ ^ qn,tQl+(nmod N),t T C ^ ^ Qn,t(^Qn,t+ 1 ~P T^ n ,t+ l)T 


n=1 


n= 1 


a N b N 

+ 2 + 2 X^ n ’ t+1 + m n,t+ 1) 2 ~ m n,t+lQn,t+l + V( q n ,t ), 

n= 1 n=1 

where = {q n ,t}n= 1 , with q n j being the coordinate of n-tli particle n = 1...N 
at the moment of time t = Z, and rri q n 4+1 , m p n t+ \ are integer numbers which stand 
for winding numbers along the q and p directions of the 2N-torus. The coefficients 
a,b,c,d are constants which will be specified below and V(q) is a smooth periodic 
function V(q + 1) = V(q). In the following we refer to the map generated by the 


action (3.1) as non-perturbed (resp. perturbed) coupled cat map in the case when 
V(q) = 0 (resp. V(q) ^ 0). 

The map shares the most important properties with the model (2.1). In 


particular, it is defined in such a way that all particles are coupled in a uniform 
way with its nearest neighbors. Models of this type are known as coupled lattice 
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maps. Their dynamical properties have been extensively studied in the last three 
decades |2Tj. Of special interest for us, is the result from |22| . where existence of 2D 
symbolic dynamics was demonstrated for a particular model of coupled lattice map. 
In the following we show that 2D symbolic dynamics with the properties described 
in Section 2.3.1 can be constructed also for the map & N . 

3.1. Dynamics. The equation of motion is generated using p n . t = —dS/dq n> t, 
Pn,t+ i = dS/dq n , t+ i: 

(3.2) qn,t- -|-1 — ( Pn,t 0/qn,t ) ~ (<?n+l,t T Qn—1,1 ) ^'n,t+l ~^ 


Pn,t +1 ( bp n j T (c (*?n+l,l T qn—l,t) 


m. 


V 

n,t+1 




Under the condition that V(q) = 0 these equations can also be written in the matrix 
form: 

Z t +1 = ® N Z t mod 1, Z t = (, qi, t ,Pi,t, ■ ■ ■ qN,t,p N ,t) T , 
with 2N x 2 N matrix given by: 


(3.3) 


23 w = 


A = — 

c 


ab 


B 

A 

B 

... 0 

0 

0 

B 

A 

... 0 

0 

0 

0 

0 

... A 

B 

\B 

0 

0 

... B 

a) 

a 

lA 



1 ( 

-c 2 

b) 

5 

B = - 

-c 


d 0 
db 0 


For the map to be continues the elements of the matrix fBjv must be integers. To 
satisfy this requirement we set c = —1 and a, 6, d G Z. 

By excluding the momentum from the system of eqs. (3.2) the dynamical equations 
for the time evolution can be cast into the Newtonian form: 

(3.4) c(q n f+ 1 + q n ,t- 1 ) + d(q n+ ij + q n -i,t) — — (o + i>)q n ,t + ™ n ,t — V'(Qn,t), 

i g n,t 


m n ,t = ~bm q nd - cm q n t+1 


+ "Cn 


3.2. Spectrum of 23^. In order to establish the type of dynamics generated by 
we need to determine the spectrum of the matrix 23 v . Since 23^ is a circular matrix 
which commutes with the shift operator its spectrum can be found explicitly. After 
substitution Z = (qi,pi ,... ,q N ,pN ) into the spectral equation 

AZ = <3 n Z 

we obtain the following condition 

IQn = Qn+l + Qn-l, cfy = -5 + A + j, S = a + 6, 
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which is equivalent to the iterative equation: 
(3.5) 


Qn +1 
Qn 


-1 

0 


Qn 

q n ~ i 


In order to close this equation the eigenvalue of the matrix above must be of the form 
e 2 mk/N^ /j — i 5 ... jv. This yields to the following quadratic equation for eigenvalues 
of *B n : 

A + A -1 = s + 2dcos(27ik/N), 


k — 1, 


N. 


(3.6) 

Accordingly, is fully hyperbolic iff |s| > 2\d\ + 2. In this case all solutions A ±(k) 
of this equation are paired such that A +(k) = Al 1 (k) and |A+(fc)| > 1 for all k. In 
what follows we will assume that this condition is always satisfied. 

3.3. Periodic orbits. In this paper we primarily focus on the periodic orbits of 
<& N . In the case of fully hyperbolic dynamics the total number A f(T) of MPO with 
period T can be easily obtained by the following formula: 


N 


(3.7) 


V(T)=|det(/-»Uj«P[AR«:), 


k =1 


where A + (k) are solutions of eq. (3.6) with the largest absolute value. It is straight¬ 


forward to see that A f(T) grows exponentially both with N and T. 

The action of periodic orbits can be, in principle, obtained by summing up all the 


terms in eq. (3.1) along the trajectory. Note, however that the generating function 
S(q t , qt+i) also contains half-integer factors b{rn q n t ) 2 /2 which are in fact irrelevant, 
as far as, the dynamical equations of concern. For convenience we omit these factors 
from the definition of actions. Such omission also does not affect the semiclassical 
theory, where only a fractional part of the actions plays a role, see e.g., pans]. 
Accordingly if T is a MPO with a period T, its action is defined by: 

t n , b \ 

(3.8) F r = EE dq n ,tqi+(nmod N),t + Cq n , t q n ,t+ 1 + 2 1 n,t + 2 Qn,t+1 + V(Qn,t) J + 

t=1 n=1 A / 

T N 

+EE ( c qn,t m n ,t+l bqn,t+l m n ,t+l ) 


f=l n =1 


where q n j, rn~ n q f\ are coordinates and winding numbers along T. Using eq. (3.4) this 
expression can be further simplified leading to: 

„ T N 


(3.9) 


Sr = -- ^ m nfi q n ,t + q n ,tV{q n ,t) ~ 2 V(q n ,t) 


t= 1 n= 1 


Remark 3.1. It worth mentioning that in most of the works on cat maps the actions 
of periodic orbits are defined in a somewhat different way, see e.g., [16]. Namely, if 
a periodic orbit T has a period T one looks for a fixed point (q, p) of the map 
and then defines t he a ction of Y as S'r(q) = <S'(q,q), where the generating function 
S is given by eq. (3.1) with the parameters a,b,c,d defined by the map rather 
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then by 4V- It can be shown that the difference between two actions Sr and Sr is 
always an integer. Since in quantum mechanics only fractional part of actions are 
of a relevance one can use both definition equivalently for semiclassical calculations. 


We, however believe that (3.9) is a more natural one from a classical point of view as 
its properties resemble ones of periodic orbit actions in generic Hamiltonian systems. 
Note, that in contrast to Sr, Sr(q) depends on the choice of the initial point q at V 
i.e., 5r(q) f Sr($ N ■ q) and it grows exponentially with T. 

4. Time-particle duality 


Remarkably, for c = d the equation (3.4) becomes symmetric under the exchange 
of times and particle numbers: n CA t. This immediately leads to the following 
duality relationship between periodic orbits of the system. 


Lemma 4.1. Let^k be a sequence of maps (3.1) with some fixed parameters a, b, c = 
d and varied number of particles k — 1,2,.... For each MPO V of & N with a period 
T there exists corresponding MPO R of the map 4> T with the period N, such that 
both trajectories run through the same set of points in the configuration space. Fur¬ 
thermore, the actions of two orbits coincide Sr = Sr' ■ 


Proof: The proof is straightforward. Given a MPO T the set of points {q n> t, rn n An = 

1.. .. N]t — 1,... T} traversed by V in the configuration space must satisfy eq. (3.4). 
By defining q ri:t = q[ n , m nit = rn' t n one can easily see that the new set {q' n t , m' nt \n = 

1 .. . ,T]t = 1,... N} satisfies eq. (3.4), as well. It remains to notice that the set 

{ q' nt ,m ' nt } uniquely defines the trajectory T' of the period N for the map 4> T . The 
equality between actions then follows immediately from eq. (3.8). □ 


Remarks 4.2. 1) The above lemma holds also for the case of perturbed cat maps 
when V(q) yl 0. 2) Note that although the two trajectories T, T 7 traverse the same 
points of the configuration space, they do not have the same set of momenta. 2) For 
N = T the set of periodic trajectories becomes self-dual. This means that each MPO 
T either has a pair MPO T 7 satisfying q nt — q' tn or satisfies the self-dual constraint 
q n ,t = qt,n f° r each t and n. The first option would actually imply that T and T 7 
traverse through exactly the same points of the configuration space. We have not 
observed such pairs in our numerical simulations. 


An immediate consequence of this lemma is the following connection between num¬ 
bers of periodic orbits. 


Lemma 4.3. The number of MPOs with a period T of the map 4V is the same as 
the number of MPOs with the period N of the map <f) T . Equivalently: 


(4.1) 


(-1)* det(J - Kf) = (-1) T det(J - 03") 


N\ 


Proof: Straightforwardly follows from Lemma 4.1 and eq. (3.7). 


□ 


Since the spectrum of matrix 03 N is known explicitly, eq. (4.1) allows to extract a 
couple of interesting mathematical identities. 
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Corollary 4.4. Let cosh (3k = s/2 — cos(2-Kk/N), cosh (3k = s/2 — cos(27rfc/T), with 
T, N being integers and s > 4 then: 


(4.2) 


N 


J^[ 4 sinh 2 


fc=i 


m 


J J 4 sinh 2 


m= 1 



After taking the limit s —> 0 this equation leads to: 


( 4 . 3 ) 

fc=l k ' m =1 k ' 

In particular, for T — 1 we obtain the well known identity: 



N 

(- 2 )*- 1 ' 


4.1. Duality between partner orbits. The duality relationship has very impor¬ 
tant implication on the correlation mechanism between periodic orbits of the sys¬ 
tem. In particular, it implies that pairs of partner orbits of a period T with the 
encounter structure given by winding numbers (wi,^) for the iV-particle system 
are in one-to-one correspondence with the pairs of partner orbits of the period N 
with the encounter structure (u> 2 , u>i) for the T-particle system. For a small number 
of particles N and large times T it is natural to expect that most of the partner 
orbits poses (0,1) type of encounters. The duality relationship shows that the op¬ 
posite regime when T is large and N is small is dominated by the partner orbits 
with inverse winding numbers (0,1). In a sense these two regimes are dual to each 
other. For instance, the number of partner orbits and their action differences must 
coincide. Although the exact duality relationship i.e., Lemma 4T holds only under 
specific condition c = d, we believe that in its weaker sense, as connection between 
correlation mechanism in the above regimes it makes sense for general parameters 
a, b, c, d. 


5. Many particle partner orbits 

The main goal of the present section is an explicit construction of partner MPO’s 
with close actions. Here we restrict our attention to the case of the non-perturbed 
maps 4 *tv with V = 0. A remarkable property of the non-perturbed cat maps is that 
their periodic orbits traverse the points of the phase space with rational coordinates 
and can be found exactly rather then as an approximation. Furthermore, their 
actions are given by rational numbers and can be explicitly evaluated, as well. 

5.1. Symbolic dynamics. As it was explained in the section ([2]), a convenient 
representation of partner MPO’s can be achieved by means of 2D symbolic dynamics 
with a finite alphabet. To introduce the symbolic representation we use as the 
alphabet the set of winding numbers rri p n t , m q n f entering the dynamical equations of 
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motions (3.2). For each moment of time we define the vector 


m n ,t = 


_ I m n,t 


m: 


n,t 


M f = (mi |t m 2 ,t ••• m N ,t\ 


Note that the whole set of these vectors from t — 1 up to t — T uniquely dehnes the 
corresponding periodic orbits. Indeed, the global winding number vector Air of a 
periodic trajectory Y can be restored from Ah by: 


(5.1) 


T 

M v = Y J ^ T N tM t. 

t= 1 


This vector can be used then in order to restore the initial conditions for the MPO: 
(5.2) Z 0 = (<B T N -I)- 1 M r . 


As a result the following matrix of winding numbers, 




/ mi,i 

m 2 ,i • 

• mjvpX 

(5.3) 

M r = 

nii, 2 

m 2,2 • 

• m V,2 



\mi,T 

m 2 ,T • 



can be seen as the unique symbolic representation of T. 

Importantly, because of the local nature of the interactions the integer vectors m n ,t 
can take only a finite number of possible values. Accordingly we obtain 2D-symbolic 
representation of periodic orbits with a small alphabet, whose size is independent 
of N. 

5.2. Construction of partner MPOs . Our construction of partner MPO’s with 
the prescribed encounter structure can be divided into few steps which are described 
below and schematically summarized in the fig. [8j 



FIGURE 8 . 4-Step construction of partner MPOs. 

Step 1 (fig. Qi). As the first step, we generate a random periodic trajectory To 
and its symbolic representation Mp 0 by means of the following procedure. Let Ado 
be an arbitrary vector of integers. Then the initial conditions for the corresponding 
MPO can be obtained by 

(5.4) 


'.M,, modi, 
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where Zq is the 2N dimensional vector of coordinates and momenta at the initial 
moment of time. The rest part of the orbit is obtained by iterative action of the 
matrix 23 at on the vector Z 0 . Note, that since the entries of the matrix 23^ are 
integers, the coordinates of the MPO are rational numbers. This property of the 
map allows to calculate the orbits explicitly avoiding any numerical approximation. 

Step 2 (fig. \3p). Recall that a partner orbit for any MPO T might exist only if it 
has (at least) one encounter. In other words certain pattern of symbols in symbolic 
representation of T must reaper for a number of times. Naturally, for small N and 
T a generic MPO does not satisfy this property and the purpose of this step is to 
prepare such an orbit. For that we choose some region of a given shape in the original 
Mr 0 , which is then copied and pasted into another location of Mr 0 . The resulting 
sequence of symbols M ro has the correct encounter structure, but, in general, does 
not correspond to any real MPO. 

Step 3 (fig. [sjcj. In order to generate 2D sequen ce co rresponding to a valid periodic 
orbit we now use Mr 0 as an initial data for eq. (5.1). More specifically, let M t := 
Mr 0 (t) be the f-th column of the matrix Mr 0 . We define the integer vector: 


(5.5) 


Mo = 53® 


T—t 

N 


M* 


t =i 


which is then used as the input for the right hand side of eq. (5.4) in order to 


generate the corresponding MPO. Denote this periodic orbit as T and its 2D symbolic 
representation as Mr- The crucial observation is that Mr and Mr 0 differ only locally 
at the places of encounter boundaries. As a result, we obtain the MPO T with the 
desired encounter structure. 

Step 4 (fig. At this stage we can construct the partner orbit T of T by rear¬ 
ranging the symbols Mp outside of the encounter regions, into a new 2D sequence 
Mf as has been explained in Section [2] Note that this new sequence of symbols is 
in general symbolic represen tation o f a real MPO T which can be straightforwardly 
recovered from Mf by eqs. ( 5.l[5.2 ). By construction f and T are partner MPO’s 
traversing approximately the same points of the phase space. 


Remark 5.1. In principle it might happen that the 2D symbolic sequence obtained 
from M r by exchange of the sequences A and B (see fig. [$[) differs from 2D symbolic 
sequence Mf of the partner orbit T by a small number of symbols. This should 
happen whenever one of the two points r t , n , T t',n' closely approaching each other in 
the phase space turn out to be encoded by two locally different symbolic subsequences. 
We however have not observed such a phenomena in our numerical constructions of 
partner MPO. think again 

Remark 5.2. One can expect that the 2D symbolic sequence obtained from Mr by 
exchange of the sequences A and B (see fig. [5]) sometimes differs from the 2D sym¬ 
bolic sequence Mf of the partner orbit T by a small number of symbols. It might 
happen, for instance, if some two points r f)n and r t / jn / are closely approaching each 
other in the phase space, but encoded by two different symbols. In our numerical 
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calculations performed for the non-perturbed maps, see Appendix A, we have not 
found evidence of such phenomena. 


Example. In the Appendix A we illustrate the above four step procedure by par¬ 
ticular calculation of the partner orbits for the case of non-perturbed coupled cat 
maps. 

5.3. Action differences. Keeping in mind applications to quantum mechanics it is 
of great importance to evaluate the action differences AS = Sr~Sf between partner 
orbits T, T. Indeed, in a semiclassical theory AS provide the weights with which 
periodic orbits of a given encounter structure contribute to spectral correlations [5]. 
It is known that in the case of single-particle chaotic systems the differences between 



c) 




(q w (i).p w (D) 




(q M (2),p/2)) 


(q n , 1 (2)>Pn.t(2)) 


q„ 


FIGURE 9. The figure (b) in the middle schematically shows two partner 
MPO together with the encounter domain and a contour inside. On the 
left figure (a) is shown the magnified caricature of the contour dC where 
the difference between actions is calculated. The right figure (c) shows the 
region in the single-particle phase space whose area contributes to eq. (5.10) 


actions of partner orbits accumulate at the encounters, where the distance between 
points is maximal. Below we show that a similar result holds for many particle 
systems, as well. 


By setting in eq. (3.9) V = 0 we obtain for the total action of T 
(5.6) 


* = 5 £ 


Qn,t'TTln,t- 


(n,t)ez& T 

Let T, f be two different MPO such that {q n ,t,Jf}n,t \( n A) £ ^ T }, {q n ,t, rhn,t\(n, t) G 


Z^ T } are the corresponding solutions of eq. (3.4). Let us calculate the difference 


between their actions accumulated in a certain region C of Z 


2 

NT* 


(5.7) 




4 c, = 2 A 

( n,t)£C 




(■ n,t)£C 

We will assume that inside of the domain C the two MPO are close to each other and 
the corresponding winding numbers coincide i.e., m n j. = for (n,t) G C. Denote 
by dC the boundary of C. Note that if C is a multiconnected domain then DC is 


composed of a number of disjoint components. Mult iplying both sides of eq. (3.4) 
for T by the coordinates q njt , and similarly eq. (3.4) for f by the coordinates q t)U 
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we obtain two equations containing bilinear forms in q n t) q nt . After subtraction of 
these equations we obtain: 


(5.8) 2(S r iSp ) ^ Qn,tQn,t+ 1 Qn,tQn,t-\- 1 A ^ ^ Qn,tQn,t— 1 Qn,tQn,t—l 

(n,i)s<9C-f (n,l)s5Cj_ 

A ^ ^ Qn,tQn+l,t Qn,tQn+l,t ~t~ ^ ^ Qn,tQn—l,t Qn,tQn—l,ti 


( n,t)£dC- 


(n,t)£dC+. 


where <9C)-, dC± and <9CA., <9CV- are pieces of the boundary going in the vertical and 
horizontal directions respectively, see fig. #)• Note that the points (n, t ) at the 
corners of the domain C belong to the both sets simultaneously. Remarkably the 
right hand side of the above equation contains only the sum over the boundaries of 
C. 

Let us now evaluate the difference between total actions of T and T. For simplicity 
of exposition we will assume that there is only one encounter of the order two. In 
such a case the whole set can be divided into three regions A, B,C separated by 
the boundary DC = dA U dB which passes inside of encounter regions, see fig. ©>)■ 
Accordingly the action difference can be written as a sum of three therms 


(5.9) 


AS = S ( r A) - + S', 


(A 


?(B) 


- 4 fl) 1 + 4 C) 


-s 


(CO 


By eq. (5.8) the differences between actions depend only on the boundary values 
of q Uit , q n ,t at dC. As a result the above expression for AS can be rewritten in a 
compact form as the sum over the boundary dC: 


(5.10) AS= ]T ASJ1.,+ ]T Asy, 

(n,l)s9C|| ( n,t)edC± 


where dC\\ := dC\UdC\, dC± := dC U<9CA are pieces of the boundary directed in 
the vertical and horizontal directions respectively. The local action differences are 
defined here as 


(5.11) A= [<fo,,(2) - «„,,(!)] felZ) - PUT] 

+ fe,((2) - 9„,i(l)] [p“,,(2)-?“,(!)] . <• e (II. A}, 

with q n> t(k) (resp. q n ,t(k )) being the coordinate at the encounter k = 1,2 of T 
(resp. f) and p“ t (/c) (resp. p^ t (k)) are the corresponding generalized momenta, 
see fig. ([9 |d). The relationship of these momenta to the coordinates depends on the 
direction a of the boundary at the point ( n,t ): 


(5.12) 


J>L(*0 : = 1n,t+l(k) - q„,t(k ), pi,(k) := q„+i,t(k) - q„,t(k). 


It worth noticing that the expression (5.11) for AS“ f can be interpreted as a local 
symplectic area of the region formed by the four points (q n;t (k ), t (k )), t (k)), 

k = 1,2 in the phase space, see fig. (j9j^). The representation (5.11) is essentially 
an extension of the similar result for the single-particle case. There the difference 
between actions of partner orbits is determined by an analogous symplectic area 
evaluated at an arbitrary point of their encounter [6]. 
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6. Symmetries 


So far we have ignored the possibility that the map <3>jv might posses a discrete 
symmetry. Such symmetries are of crucial importance for the associated quantum 
problem, where they determine which universality class system belongs to. On 
the classical side of the problem the presence of symmetries gives rise to partner 
orbits with close actions which do not traverse approximately the same points of 
phase space, but rather two sets of points related by the symmetry operation. For 
instance, in the systems with the time reversal invariance there exist also partner 
orbits which traverse (approximately) the same points of the configuration space but 
have opposite momenta at some pieces of the trajectories. For this to happen it is 
sufficient to have just one encounter, which is the case of the original Sieber-Richter 
pairs [7j. 

All symmetries of Hamiltonian systems can be divided into two classes. The ones 
which preserve the form of Hamiltonian’s equations are called canonical. For 2N 
dimensional cat maps they are represented by 2N x 2N integer matrices £, Z —» £Z, 
Z = {qi,Pi,q 2 ,P 2 , ■ ■ ■ ,qN,PN) J satisfying 


( 6 . 1 ) 


£03 vC 1 = 03 


N, 


<r = i. 


For instance, all maps (3.2) commute with the matrix £ = —I which effects the 
transformation q—» — q, p—> — p. On the other hand a Hamiltonian system might 
also have anticanonical symmetries, which reverse the signs of Hamiltonian’s equa¬ 
tions and the Poisson brackets. The most prominent example of the last one is time 
reversal symmetry. The anticanonical symmetries of $ are represented by 2N x 2N 
integer matrices £, that satisfy 


( 6 , 2 ) 


£»„£-' = <8 


N ’ 


£ 2 = I 


det £ = —1. 


Its easy to check that for b = 0 the map $ has the anticanonical symmetry: 


£(gi, Pi, 92 , P2, ■ ■ •, Qn,Pn) j = (pi, qi,P2, 92 , • ■ ■, Pn, 9jv) t , 


which amounts to exchanging of the corresponding moments and coordinates. 

Assuming that the map <£># posses a symmetry Sc (either canonical or anticanoni¬ 
cal) defined by a 2N x 2N matrix £ with integer entries let us describe the structure 
of partner orbits with one encounter. First observe that each MPO T has a pair 
orbit r* = Sc ■ T which is obtained by action of the linear map (£ on the set of points 
T. Both orbits have exactly the same actions. The corresponding partner orbits are 
obtained then by a recombination of different patches from T and T*. 

To be more specific consider an orbit T having one non-trivial encounter with zero 
winding numbers (0,0). Such an encounter corresponds to two frame-like regions 
Ei, E 2 of Z 2 t where 2D symbolic representations of F and its conjugate T* are 
related by the symmetry operation Sc in the following way: 


(6.3) E :— ip r (Ei) — ^r*(E 2 ), E* :— ^ r (E 2 ) — ^r*(E\). 

The 2D symbolic representation Mf of the partner orbit f is then obtained by 
copying the interiors A, B of the encounter regions from Mr* and pasting them 
into the corresponding places of Mr, as shown in fig. 10, The partner orbit of T* 
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FIGURE 10. On the two left figures are schematically shown the 2D sym¬ 
bolic representations Mr, Mp* of some MPO T and its symmetric counter¬ 
part r* having one encounter with (0,0) winding numbers. Such conjugate 
pairs of MPO’s with exactly the same action exist whenever the system 
posses a symmetry. On the two right pictures are depicted the symbolic 
representations Mp Mf* of the corresponding partner orbits F, F*. They 
are obtained from T and T* by copy pasting of the interior regions of en¬ 
counters. In general F, F* (resp. T, T*) have close but non-identical actions. 
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FIGURE 11 . On the two left figures are schematically shown the 2D 
symbolic representations Mp, Mp* of symmetrically conjugate MPO’s with 
one encounter having (1,0) winding numbers. On the two right pictures 
are depicted the symbolic representations Mp Mp* of the corresponding 
partner orbits F, F*. They are obtained from T and T* by copy pasting of 
the corresponding interior regions of Mp, Mp*, as shown in the figure. The 
partner orbits for MPO’s with encounters of the type (0,1) are constructed 
in the same way by exchanging particle and time directions. 


can readily be obtained by applying symmetry operation to F i.e., F* = ■ F. In 

Appendix B we illustrate the above procedure with an explicit construction of a 
quadrapole of periodic orbits T, T*, F, F*. 

For MPO’s with one encounter having winding numbers (1, 0) and (0,1) the con¬ 
struction is analogous. In this case the corresponding encounter regions of are 
two strips winding torus in the particle (resp. time) direction and separating it into 
two pieces. The symbolic representation of partner orbits is obtained from Mp, Mp* 
by copy pasting of the symbols from these regions, as shown in fig. 11 Note that 


for the encounters of the type (0,1) the obtained partner MPO’s are the original 
Sieber-Richter pairs in 2N dimensional phase space of the map $at. 
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7. Partner MPOs for perturbed map 


It is well known that the spectra of quantum cat maps essentially differ from ones 
of generic chaotic systems [IB]. From the semiclassical point of view this anomaly 
can be attributed to the large degeneracies among the actions of the periodic orbits 
which are given by rational numbers for a cat map. These degeneracies grow expo¬ 
nentially with the period of orbits ra leading to the distortion in the distribution of 
action differences. The genericity, however, may be recovered by adding a periodic 
nonlinear perturbation [18]. From this prospective it is natural to inquire whether 
the correlation mechanism between the partner MPOs described in the previous sec¬ 
tions carries over to the perturbed coupled cat maps. Below we demonstrate that 
the answer to this question is indeed affirmative. 


Consider the model of perturbed coupled cat maps generated by the function (3.1), 
where some periodic potential V(q) = ft Vo (q) is turned on. We require that, 
maxq |Vo(g)| = 1 and relate to ft as the strength of the perturbation. The par¬ 
ticular choice of the parameter ft is dictated by the requirement that the perturbed 
map <Fjv remains fully chaotic. By the structural stability of classical chaos this is 
guaranteed for sufficiently small ft. A maximum possible range for ft can be esti¬ 
mated by considering the linearized form of the map Let q n> t,n = 1,2,... ,N 


be some solution of eq. (3.4) and let 8q n j be a small deviation from it. After sub¬ 
stituting q n ^ + 5q n<t into the equations of motion (3.4) and leaving out only linear 
terms we obtain 


(7.1) 8qn,t- i-i T 8q n ^t—\ d(5q n ^.\^ -\- 8qn—i,t) (s V 

By comparison with the non-perturbed case analysis in Section 3.2 we arrive to the 
inequality \k\ < fc max = |s| — 2|rf| — 2. It can be expected that under this condition 
the perturbed map <& N is fully hyperbolic. 


7.1. Construction of MPO partners. In general, calculating periodic orbits of 
a non-integrablc system is a non-trivial task. To this end a number of methods 
have been developed, see e.g., [24]. Our approach is based on the minimization 
of the perturbed action, i.e., the generating function S(q n j) with respect to q Uit . 
The procedure goes as follows. One starts from a MPO r 0 = {q^j \ (n, t) e Z^ T } 
of a non-perturbed map and adds a small perturbation 8nVo(q) to the action. By 
structural stability of chaotic dynamics it can be expected that, for a small 5k the 
non-perturbed orbit lays in the vicinity of the perturbed one r x = {q^} \ (n, t ) G Z^ T } 
which can be easily found by minimization of the action S(q n j) by using the standard 
gradient descent method. At the next step one takes {//„,} as the initial data for 
the gradient descent method and obtains a new MPO T 2 for the doubled strength of 
the perturbation 25k. By iterating this procedure m times one arrives to a periodic 
orbit Y m which is a smooth deformation of T 0 for the perturbation of the strength 
k = mSK < K max . Importantly, if we start from two partner orbits r 0 ,r 0 with the 
related matrices M ro ,Mf 0 of winding numbers, this relationship persists at each step 
of the iterations. This means that the resulting MPOs T m , Y m are again partners 
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with small action differences, whose symbolic representations are related in the same 
way as ones of the initial orbits. 

In Appendix C we give an example of explicit calculation of a pair of MPOs 
T, r for the perturbed coupled cat map and demonstrate stability of the pairing 
phenomena. Namely, we show that all the points of T and f are paired and their 
action differences accumulate in the encounter regions. 


7.2. Action differences. It is instructive to see how the formula (5.10) for action 


differences between the partner orbits is affected by the perturbation. In accordance 


with eq. (3.9) the actions get extra terms depending on V. After repeating the same 


steps as in the derivation of the non-perturbed case we end up with the additional 


term in the right hand side of eq. (5.8): 


§s {c) = ^2 A f V(qn,t) + V(q n , t )) ( q n ,t ~ q n ,t) ~ V(g„, t ) + V(g n , t ). 

(n,t)eC 1 ' ' 


Like the main contribution in eq. (5.8), the term dS^" 1 is accumulated mainly in the 
encounter region. It is straightforward to see, however, that 5S <yC ' 1 is proportional 
to the third order of 5q n . t := q n , t - q Ujt : 


6S^ C] = ^ 


5q 3 


n.t 


12 


V {3) 


qn,t T qn,t 


+ 0(Sql t ), 


while the principle term S'p 6 ' 1 — S^" 1 is of the order Sq^ t . As a result, SS^ can 
be dropped out for MPOs with large encounters (where the distances between the 


partners are small). Since the same conclusion holds for all other terms in eq. (5.9) 


it follows immediately that the formula (5.10) is still valid in an approximate sense. 


The larger the size of encounters, the better approximation gives the equation (5.10). 


8. Conclusions 

8.1. Quantum problem. So far we have studied clustering of MPOs in the con¬ 
text of purely classical theory. While leaving applications to quantum mechanics 
for future investigations, let us briefly discuss the implications of our results for 
a semiclassical theory of many-particle systems. Spectral correlations in quantum 
systems with fully chaotic classical dynamics can be related to correlations between 
periodic orbits. Consider as an example the spectral form factor K 2 which is defined 
as the Fourier transform of the two-point spectral correlation function. Via the trace 
formula it can be expressed as a double sum over pairs of partner periodic orbits. 
These pairs r,f can be clustered in accordance to their topological structures Mr, 
such that: 

(8.1) K 2 (T,N) = Re E E M r exp(iASr/h ef f), 

t { r,r|j- r =j-} 

where the first sum runs over all possible topological structures T of MPOs [25] . 
The factors Mr, ASr are determined by stabilities and by differences between actions 
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of partner orbits, respectively, with being the integer Hilbert space dimension 
of the corresponding quantum cat map, see [lijj . 

In order to clarify which topological structures of MPOs are of a relevance for 
given number of particles N and time T, it is of crucial importance to estimate 
how many periodic orbits T posses a given topological structure T. The amount of 
such orbits is determined by the number of possibilities to choose the encounter sets 
En\ k — 1,.. .l n within the prescribed topological order. As we argue below this, in 


Note that 


turn, depends on the winding numbers w n of En^s and their widths p r 
systematic contribution into the sum (8.1) comes from the terms for which AS'r is 
of the same order as h e ff. This fixes the scale of encounter width as p n ~ | log h e ff |. 
Consider now an encounter set having the winding numbers w = (0,1), i.e., winding 
around the torus in the particle direction. Such a set occupies at least N ■ p points 
of Z 2 T . Analogously, for encounters with winding numbers w = (1,0), w = (0,0) 
the minimal number of points is given by T ■ p and p 2 , respectively. It is easy 
to see that the number of periodic orbits with the fixed symbols at the encounter 
domains En \ k = 1,... l n grows exponentially with the total number of points of 
Z^ T outside of the encounter sets. From this we can immediately conclude that 
for N -C \ \og h e ff \ <C T most of the periodic orbits have encounters of the type 
w = (0,1), since they occupy the minimal “area” of Z^ T . For such parameters N, T 
the correlations between actions of MPOs can be understood just by interpreting 
the system dynamics as single-particle motion in many-dimensional phase space. In 
other words in this regime the semiclassical theory of many-particle systems should 
be the same as for single-particle systems with chaotic dynamics. In the regime T <C 
| \ogh e ff \ -C N the dominating encounter type of periodic orbits is w — (1,0) and 
the clustering mechanism is dual to the previous one. This means that correlations 
between actions of periodic orbits can be accounted by utilizing single-particle theory 
with the parameters N and T exchanged. Finally in the case when both parameters 
are “large”: T,N | log^ e //1 most of the periodic orbits posses encounters with 
zero winding numbers w = (0,0). In this case the correlations between actions 
of MPOs cannot be accounted on the basis of single-particle interpretation of the 
system dynamics. For this regime a new genuine many-particle semiclassical theory 
has to be constructed. 

8.2. Generalizations. The theory proposed in this paper has been constructed 
and verified in the framework of a specific model of coupled lattice maps. One 
might wonder whether it carries over to the Hamiltonian systems with continues 
time evolution. We believe that the main results hold also for generic Hamiltonians 
of the form (2.1) provided that the dynamics are chaotic. In order to clarify this 


point consider, as an example, a many-particle periodic orbit T(t) = {{x n (t),p n (t))} 
where t £ [0, T] now continues time and the coordinates x n (t) satisfy the system 
of N Newtonian equations generated by the Hamiltonian (2.1). Assume that T(£) 
has an encounter domain with the frame-like structure shown in fig. [7} In the 
coordinate form this means that for (n, t) £ E we have an approximate equality 
x n (t) ~ x n+n '(t + tf), where (n + n',t + t') £ E' are the points of the shifted 
encounter set. (For sufficiently large T and N such encounters should exist just by 
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statistical reasons.) Let us now construct another trajectory following the rules for 
search of partner orbits. To this end define: 

(x n (t) if (ra, t) (£ B |J A, 

x n (t) \= lx n+n \t + t') if (n,t)eB, 

[x n - n >(t-t') if (n,t) e A. 

It is now a straightforward observation following from the locality of the Hamilton¬ 
ian interactions that the set of trajectories T(t) = {(x n (t),p n (t))} approximately 
satisfies the equations of motion. It follows then from the shadowing lemma for 
flows of chaotic systems ra that a real periodic orbit close to T(f) should exist. By 
construction this new orbit is the partner of T(t). 

Taking a step further one can also question whether discretization in the particle 
direction might be left out as well. Making n continues naturally leads to a field 
model, where the held 0(x, t ) satisfies some homogeneous differential equation and 
the continues parameter x G [0, L\ plays the role of n. Assuming that the resulting 
dynamics are chaotic it can be expected that periodic solutions cp per (x,t) of the 
differential equation posses similar properties as MPOs in the coupled lattice maps. 
In particular we expect that the mechanism of correlations between MPOs carries 
over to (j) per (x,t ). We believe that this correlation mechanism can be of importance 
for development of the semiclassical approach to the corresponding Quantum Field 
Theory. 

Another natural extension of the present setup arises if the particles are ordered 
with respect to the d-dimensional lattice i.e., n G Z d , d > 1 rather then to the one- 
dimensional Z. In this case propagating MPOs sweep d + 1-dimensional “surfaces” 
in the one-particle phase space. It would be of interest to study how the dimension 
d affects the correlation mechanism between MPOs. 

8.3. Summary. Periodic orbits of chaotic systems can be organized into clusters 
of partner orbits where all elements have approximately the same actions. The 
clustering phenomena plays a crucial role for the semiclassical theory. So far how¬ 
ever, it has been studied only on the level of effectively “single-particle” systems. 
In the present work we considered clustering mechanism among periodic orbits of 
fully chaotic many-particle systems which dynamics are governed by translation- 
invariant Hamiltonians with local interactions. On the heuristic level periodic orbits 
in such systems can be visualized by discretized two-dimensional “surfaces” in the 
one-particle configuration space. A more rigorous description is achieved by using 
symbolic dynamics. Contrary to the single-particle formalism, where periodic or¬ 
bits are encoded by linear symbolic sequences, here the symbolic representation is 
provided by 2D toric arrays of symbols. A periodic orbit T has a partner f when 
its symbolic representation Ar has at least one /-encounter - a 2D subsequence of 
symbols repeating itself for l > 1 of times. The symbolic representations Af of the 
partner orbits are then obtained by rearrangement of symbols Ar in accordance to 
certain prescription. The differences between actions of T and its partners accumu¬ 
late mainly in the encounter regions and decrease exponentially with their width. 
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Different families of partner orbits can be distinguished by the corresponding topo¬ 
logical structure diagrams. Each structure is represented by a toroidal surface glued 
to itself along a number of one-dimensional lines representing the encounters. The 
diagrams with encounters winding around torus torus only in the particle direction 
are the ones which obtained in the single-particle interpretation of classical dynam¬ 
ics. The other structural diagrams are essentially different. They represent partner 
periodic orbits whose correlation mechanism differs from one found in single-particle 
systems. These new families of partner orbits seem to be essential for construction of 
a proper semiclassical theory in the case when the number of particles is sufficiently 
large. 
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Figure 12. Symbolic representations Ap (left) Af (right) of the two 
partner periodic orbits (r, f) for parameters a = 3, b = 2, c = d = —1, 
T = 50, N = 70. The corresponding elements of the matrices Mp and 
Mp are encoded by the hexadecimal numbers in accordance with the 
table ( |8.2 ). The identical encounter regions are highlighted by green 
(color on-line). The coinciding inner regions of the encounters (which 
are exchanged in T and T) are colored by the identical colors (blue 
and dark blue). 


Appendix A 

In this section we construct explicit example of two partner MPOs for the map 
& N with the parameters: N = 70, a = 3, b = 2, c = d = —1. The initial vector 


M 0 in the equation (5.4) is generated as a set of random integers taken from the 
interval [—10 21 ,10 21 ]. The two partner MPOs of the period T = 50 were generated 
following the protocol described in Section 5.2, see also the explanatory scheme 
on the figure (§]). The resulting symbolic representations Ap, Ap are presented on 
the picture 


12 


where each of 16 valid pairs of winding numbers m = ( m,^ q \m ^) 
see figure 12h is encoded by a symbol a from the hexadecimal alphabet A. The 


correspondence between symbols of A and winding numbers is provided by the table 
below: 


( 8 , 2 ) 


lU 



mmm 



3 4 


6 7 


9 A B CD 
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FIGURE 13. Left: The phase-space coordinates of the pair of MPOs T, f 
encoded by the symbolic 2D sequences Ar, Af on the fig. 12 The positions 
of particles belonging to T are depicted by the green circles, while the 
positions of particles of the orbit T are the blue circles of smaller sizes. 
All green and blue points come in pairs and there are no unpaired points. 
Right: The enlarged region of the phase space (the selected rectangle on 
the left picture). The geometrical centers of the circles visually coincide 
for almost all points except the four points in the middle of the diagram. 


On the fig. 12 the encounter regions of the orbits Ar, Af are colored in green. One 
can easily verify that the symbolic sequences inside these regions coincide with each 
other. Note that the shape of the green areas slightly differs on its borders from 
the rectangular shape (colored in yellow) of the original encounter-patches. As it 
is clearly seen, the deformation on the border of legitimate fields causes only local 
transformation in the symbolic representation. 

To evaluate the actual distances between the partner orbits in the phase space 
we have restored their coordinates and momenta T = {{q n ,t,Pn,t)\( n , t) G Z^ T }, 
f = {( qn,t,Pn,t)\{n,t ) G Z)( T } from the corresponding matrices of winding numbers 
M r ,Mf. The resulting sets of the points in the phase space V are shown on fig. 13 
by the black points for T and by the green circles for its partner f. As one can 
immediately see, all the points of two orbits are coupled into pairs separated by 
very small distances. For some points the distances are somewhat larger and the 
pairs form quadruples, see the enlarged region of the phase space on the fig. 13 
These special points correspond to the encounter region of the orbits. 

To quantify the metric distances between partner orbits we have calculated the 
mean square displacement in the phase space for each couple of the paired points. 
The value of the displacement d n j for the particle number n at time t is given by 


( 8 . 3 ) 


d"n,t \J ( Qn,t Qn' ,t')^ H - (Pnjt Pn'i 
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Figure 14. Diagram of metric distances between MPOs F and T 
from the fig. [12} The distances are calculated according to the for¬ 
mula (8.3). The maximal distance between points is ~ 1.7889 x 1CT 3 . 
The numeric error is of the order ~ 10 -14 . 


where (<V ,t', Vri ,t') G T is the point paired to the one ( q n ,t,Pn,t ) £ T. The diagram of 
the metric distances is plotted on the figure 14 As expected, the largest distances 
d U)t ~ 10~ 3 are observed in the encounter regions (appearing as white areas on the 
figure 12). Outside of the encounters the separations between partner orbits are 
extremely small with d nt t < 1CF 12 . 


Appendix B 

In this appendix we explicitly construct the quadruple of orbits T, T*, f, f* of 
the period T = 50 for the same map 4V as in Appendix A. As has been explained 
in Section [6] this map posses the symmetry S^: q —> 1 — q mod 1 , p —y 1 — pmod 1 . 
To construct partner orbits which traverse different points of the phase space V 
we utilize 4-step protocol similar to one used in Appendix A. First, we generate 
a random MPO r 0 , its symmetric counterpart Tq and the corresponding matrices 
of winding numbers Mr 0 , Mr*. We then take a frame-like region E\ C Z 2 T and 
superimpose the local winding numbers of Mp* at E 1 upon another (shifted) region 

E‘2 C Z 2 t of Mr 0 . The resulting matrix M ro has two regions E\ , E 2 satisfying 
Note, however, that 


(6.3) 


Ir 0 does not correspond to any real MPO. In order to 


obtain matrix of winding numbers for a valid periodic orbit we then apply eq. (5.5) 
(see step 3 in Section 5.2). This yields to the matrix Mr of winding numbers for a 
valid trajectory T with the correct encounter structure. Its symmetric counterpart 
is obtained by application of the symmetry operation Sc. Finally the partner orbits 
r, f* are obtained by rearrangement of symbols from Mr, Mp* as shown in fig. 10 
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Figure 15. Symbolic representations Ar (upper left) Ap (upper 
right) Ar* (down left) Ap* (down right) of two partner periodic orbits 
(r, f) and their symmetrically conjugate counterparts (r* = Sg • T, 
r* = Sc • r) for T = 50, N = 70. The identical regions of symbols are 
highlighted by the same colors (color on line). The two encounter re¬ 
gions (green and dark-blue) are pairwise conjugated by the symmetry 
operation Sf. 


The resulting symbolic representations of orbits T, T*, T, T* are shown in fig. 15 
The elements of the matrices Mr (resp.) Mp are encoded here by the hexadecimal 
numbers in accordance with the table (8.2). In order to see that T, T*, T, f* 
are indeed partner orbits we plotted on the figure (16) the corresponding sets of 
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FIGURE 16. Left: The phase-space coordinates of the quadruple of 
orbits T, T*, T, r*. The positions of the particles of the trajectories T 
(T*) are depicted by the large dark (light) green circles. The phase space 
coordinates of the partner MPOs f (f*) are depicted by blue (red) circles 
respectively. All points come in pairs and there are no unpaired points. 

Right: The enlarged region of the phase space (the selected rectangle on 
the left picture). The geometrical centers of the circles visually coincide for 
almost all points except the four points in the middle of the plot. These 
points belong to the encounter region, where the distances between partner 
orbits are maximal. 

momenta and coordinates. It is clearly visible on this plot that any point from F 
(resp. F*) closely approaches one of the points from either T or T*. In other words 
all the points of f and T are paired up to the symmetry operation Sg. As expected, 
the largest distances between the pairs of the points are observed in the encounters. 

Appendix C 

In this appendix we present explicit calculations of partner MPOs for the per¬ 
turbed map. The periodic potential V(q) is chosen in the form 

hz 

V(g) = -^cos(2vrg), 

with the other parameters of the map set to the same values as in Appendix A. For 
the sake of convenience we require that the alphabet of the symbolic representation 
(e.g., the number of symbols) would not change under the perturbation. This can be 
achieved by a proper choice of the perturbation strength, which in our calculations 
was set to k — 1. As we found in our numerical study for this choice of k the 
dynamics of the system are chaotic and the symbolic alphabet is still preserved. 

The initial partner MPOs To, To of the non-perturbed map were taken from our 
calculations in Appendix A. The partner orbits T, F for the perturbed map were then 
constructed following the recipe in Section [TJ The two obtained sets of points are 
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FIGURE 17. Left: The phase-space diagram of the pair of MPOs in the 
case of perturbed map. The coordinates of the orbit and its pair are de¬ 
picted by the green circles and by the blue circles of smaller size, respec¬ 
tively. All green and blue points come in pairs and there are no unpaired 
points. Right: The enlarged region of the phase space (rectangle on the left 
picture). The geometrical centers of the circles visually coincide for almost 
all points except for the four points in the middle of the diagram belonging 
to the encounter. 



FIGURE 18. 3D diagrams of metric distances for non-perturbed map (left) 
(2D version of this figure is shown on the fig. 14) and for the perturbed 
map (right). The maximal distance in the case of the perturbed map is 
~ 1.90578 x 10~ 3 , the numeric precision is of the order ~ 10 -10 . 


depicted on the fig. 0 As one can see, all the points of partner orbits are coupled 
into pairs separated by very small distances. In the encounters the distances become 
larger and the points form quadrapole, see the enlarged region of the phase space 
on the fig. 17 The metric distances between T, T were evaluated using the same 
method as in Appendix A. The result is presented on the fig. 18 in 3D form. For 
a comparison, the result for the original non-perturbed MPO’s To, To is shown on 
the same figure. Although the points of the initial r 0 , f 0 and perturbed orbits T, 
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T are positioned far away from each other, the generic picture of encounters has no 
structural changes. One can see only slight changes in their pairwise differences on 


the fig. 18 
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